Breakdown of the universality in few-boson systems 

Arnoldas Deltuva 1 '!*] and Rimantas Lazauskas 2 '^ 

1 Centro de Fisica Nuclear da Universidade de Lisboa, P-16^9-003 Lisboa, Portugal 
2 IPHC, IN2P3-CNRS/Umversite Louis Pasteur BP 28, 
F-67037 Strasbourg Cedex 2, France 
(Dated: September 7, 2010) 

Abstract 

We develop a series of resonant short-range two-boson potentials reproducing the same two-body 
low-energy observables and apply them in three- and four-body calculations. We demonstrate that 
the universal behavior predicted by effective field theory may be strongly violated and analyze the 
conditions for this phenomenon. 
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I. INTRODUCTION 



Effective field theory (EFT) found a great success in particle, nuclear, and condensed 
matter physics. EFT and the alternative quantum-mechanical approach with EFT-generated 
interactions turn out to be a very practical tools to deal with the problems involving largely 
separated energy (or length) scales. One such example, revealed by the nuclear and cold 
atomic gas structures, is a system of interacting particles at low energies whose behavior 
is governed by the large two-particle scattering length a , i.e., |a | ^> r , r being the 
range of the interaction, while by low energy we understand that characteristic center-of- 
mass (cm.) momenta k are significantly smaller than l/r ). In this case, in order to 
describe the occurring physical phenomena, one naturally introduces EFT expansion in 
small dimensionless parameters ro/ao and Zero, which should be well behaved apart from the 
well-known case of singular zero-range interaction (r — > 0), leading to Thomas collapse in 
the N > 2 particle systems with attractive pairwise interactions. 

The r > case has been first seriously investigated by Efimov jl, 2], who has derived 
several interaction-independent scaling relations for three-particle observables of the systems 
with the particles constrained in the region where their interdistance R satisfies the relation 
ro <C R <C |ao|. Notably, as the ratio oo/to increases, the accumulation of the weakly 
bound three-particle states (Efimov states) has been demonstrated. Efimov states have 
found much attention recently. Properties of the simplest many-particle structures, namely 
+ states in the systems of identical bosons, have been classified. In particular, numerous 
universal (interaction independent) relations have been established between three-boson 
observables. It has been accepted that in the first approximation (leading order of EFT) all 
the low-energy properties of three-boson system are set by two parameters: one two-body 
parameter — two-boson (dimer) binding energy or scattering length, and one three-body 
parameter — three-boson (trimer) binding energy or particle-dimer scattering length {3, 4]. 
More recently, it has even been demonstrated that properties of the four-boson system 
(tetramer) are governed by the same two parameters and no additional four-body scale is 
required to establish universal relations between three- and four-boson observables {5, Q]. 
Universality studies have been extended also to fermionic three- and four-body systems, 
by means of the so-called pionless EFT, recovering phenomeno logically observed Phillips 
(triton binding energy vs. neutron-deuteron doublet scattering length) and Tjon (triton vs. 
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alpha particle binding energies) correlations [7J. Most of these relations have been derived 
relying on the purely attractive (and often contact) two-particle interactions, while N > 2 
systems are balanced to prevent Thomas collapse by introducing repulsive three-particle 
force. Such approach evidently helps to preserve Efimov relation r <^ R ^ \a \ in all 
generated structures. Nevertheless, in real physical systems, if the two-particle cluster is 
resonant, one can guarantee only R <C |ao| condition but one cannot always assure that 
three- or four-particle system prevails the condition ro <C R, i.e., most of the time particles 
stay outside the interaction region in the few-particle system. 

In this work we follow the quantum-mechanical approach with slightly different interac- 
tions, namely, we develop series of short-range separable rank 2 potentials preserving the 
same two-boson binding energy and low-energy scattering observables. These potentials 
have been derived somehow mimicking interaction between 4 He atoms, i.e., choosing the 
boson mass h 2 /m = 12.12 K-A 2 and fixing the two-boson scattering length to a =104.0 
A, but we admit that a realistic description of 4 He multimers is not the aim of the present 
work. Instead, with these potentials we calculate and compare series of three- and four- 
particle properties, including binding energies, mean square radii, and scattering lengths. 
We demonstrate that, despite resonant two-boson interaction, the breakdown of the univer- 
sality and commonly established EFT results is possible, and investigate conditions for the 
occurrence of this phenomenon. 

In Section [Til we describe the dynamic input and recall the methods used to solve few- 
particle equations. In Section UTTJ we present the results for three- and four-boson bound 
states and scattering lengths. Section [TV] contains our conclusions. 



II. DYNAMICS 



As a dynamic input to our calculations we use a set of two-boson potentials that reproduce 
the same two-particle binding energy and low-energy scattering observables. Such a set could 
be generated from the initially chosen potential v < -°- ) by unitary transformations like those 
of the renormalization group [8]. The resulting potentials are nonlocal even if the initial 
one was local. We follow a simpler procedure and use separable (and thereby nonlocal) 
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potentials 

v = Y^\9i)Kj{j3j\ (!) 

that, for simplicity, are restricted to act in the S'-wave only. In close analogy with the 
most EFT-generated potentials [5J, the form factors \gi) are chosen as Gaussians that in the 
momentum-space representation are defined as 

(klm) = e- (fc/A ' )2 (2) 

whereas in the configuration-space they are 

(r|^) = (A J /v / 2) 3 e-( A ^ 2 . (3) 

The initial potential was chosen to be of rank one, i.e., Xff = SuSjiXf'i , with the 
momentum cutoff = 0.2 A -1 and A^ = [irm(l/ao — A^/\/27r)/2] _1 determined by the 
two-body scattering length a =104.0 A. The corresponding dimer binding energy and the 
effective range are B 2 = 1.318 mK and r = 15.0 A, respectively; the latter is about two times 
larger than the value provided by the "realistic" interaction models for 4 He atoms. All the 
other potentials were obtained by choosing Aj = jA with A values in the range between 0.133 
and 0.305 A -1 , i.e., 13.8 < a A < 31.7, and determining the respective strength parameters 
Xij from the fit of the calculated observables, i.e., B 2 , ao, and two-particle scattering phase- 
shifts up to about 50 B 2 cm. energy, to the predictions of the initial potential v ( -°- ) ; we found 
that rank n r = 2 potential with three free parameters An, A22, and A12 = A21 is sufficient 
for a high quality fit with a typical four digit accuracy. Furthermore, a broad interval of 
the scattering energies included in the fit, roughly corresponding to the natural energy scale 
H 2 /mrQ [5], guarantees that in the effective range expansion not only k 2 but also higher 
order terms are well retained. 

The separable form of the potential is convenient for the fit but the advantage of the sepa- 
rability is not used when solving few-body equations that can include any form of the poten- 
tial. The exact three- and four-body bound state and scattering equations are solved numer- 
ically by two completely different methods. The symmetrized form of equations is employed 
which is appropriate for the system of identical bosons. The configuration-space Faddeev 
and Faddeev- Yakubovsky (F/FY) integro-differential equations [9|JlO( for the wave function 



components are solved using the numerical technique of Ref . [ll| . In the momentum-space 
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framework the integral form of F/FY equations is used for the description of the bound 



states whereas in the case of scattering the A 
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131 ] for the transition operators (see Ref. 



t, Grassberger, and Sandhas (AGS) equations 



14j for the symmetrized version) leading more 



directly to the observables are solved following the technical developments of Refs. |15l . 116 ] . 
Both calculations are performed in the partial wave basis. Although the two-body interac- 
tion is limited to the S wave, higher angular momentum states up to l y , l z < 2 are needed for 
the particle-pair and particle-triplet relative motion. With those states included we obtain 
well converged results, while both methods are in excellent agreement, the differences being 
well below 0.1%. 



III. RESULTS 

We would like to note first that 4 He multimers in the quantum-mechanical approach with 
EFT-generated potentials have been also studied by Platter et al. [5]. However, rank n r = 1 
potentials have been used in the latter work with An fitted to the dimer binding energy 
B2; although in the leading order in (aoA) -1 the relation ao = [^/(mf^)] 1 ' 2 holds, strictly 
speaking obtained ao as well as the scattering phase shift values are cutoff-dependent in this 
approach. To renormalize the three-body system, i.e., to get (nearly) cutoff- independent 
results for physical three-body observables, an additional separable three-boson force (re- 
pulsive or attractive, depending on the A value) had to be included in the calculations of 
Ref. j^]; otherwise, the three- and four-boson binding energies show quite strong dependence 
on A which is interesting to compare with our rank 2 results. We perform n r — 1 calcula- 
tions with An fitted to the same two-boson scattering length ao = 104.0 A as in the n r = 2 
case. The two approaches provide very different A-evolution of the few-boson binding ener- 
gies as Fig. [T] demonstrates. The ratios of trimer and tetramer binding energies B% ] to the 
dimer binding energy B 2 for the rank 1 potentials increase monotonically with A, whereas 
for rank 2 they show much weaker dependence in the range aoA < 26.0 where all £?jy fB 2 
slowly decrease with A. In that region all potentials support ground and one excited state 
for trimer and tetramer, whereas the tetramers excited state binding energy B± is only 
slightly larger than B^ and thus follows closely the trimers ground state curve; the same 
applies to the n r — 1 case in the whole considered region. From rank 2 calculations in the 
range a^A < 26.0 it may seem that when the two-boson observables are fixed by a proper fit 
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there may be no need for a three-body parameter to get (nearly) cutoff-independent results 
for the three-body and four-body system. However, for larger A values all B^ / B 2 in the 
n r = 2 case show a nontrivial behavior despite the fact that all potential parameters vary 
smoothly and slowly with A. First, at a^A = 26.5 the tetramers excited state binding energy 
detaches from B^ and exhibits a very rapid increase but stabilizes for a moment when 
it almost attains the tetramer ground state binding energy B± that, with a tiny delay, also 
starts to increase rapidly. Furthermore, at a A = 26.75 the second excited tetramer state 
appears with the binding energy B^ only slightly larger than B% . In this region of drastic 
tetramer variations the trimer properties remain practically unchanged, suggesting that not 
a three- but four-body parameter would be needed for the renormalization of the considered 
few-particle system [17J. The need of the four-body parameter as well as the existence of 
three tetramer bound states disagrees with the universality predictions {5, O]; the reason 
for this nonstandard behavior will be explained later. The trimer also shows qualitatively 
similar behavior just at larger A values. Around aoA = 27.5 the binding energies of the 
trimer ground and first excited state, B$ and B^\ respectively, start to grow very rapidly; 
the former increases much like £4 while the latter stabilizes for a while. As B^ increases 
the second excited tetramer state disappears at aoA = 28.14 and the binding energy B^ 
of the first excited tetramer state comes close to -83°^ again. Furthermore, at aoA = 28.12 
the second excited trimer state appears. Unlike the stronger bound trimers, this one has 
properties of the true Efimov state, i.e., it slides under the dimer threshold if the two-boson 
interaction is made stronger. The second excited trimer state loses that property above 
aoA = 30.7 where its binding energy B^\ much like B^ at aoA = 27.5, enters a phase 
of a rapid increase. However, the Efimov property is recovered in the third excited trimer 
state appearing at aoA = 31.7. It is interesting to note that the second excited tetramer 
state shows a similar behavior, i.e., it also disappears if the two-boson interaction is made 
stronger. Furthermore, at A values corresponding to the appearance of the second and third 
excited trimer states the binding energy of a deeper lying trimer state is almost the same in 
both cases, about 5.6B 2 . This value is slightly smaller than one found in Ref. 4| for zero- 
range potential, which should be consequence of the finite range. As can be seen from Fig. 
31 of Ref. [3]], zero-range calculations with fixed B 2 and B^ 1 slightly underestimate B^ +1 ^ 
values compared to the results obtained with realistic potentials. Therefore the critical value 
of B^ /B 2 is lower for the finite-range potentials. 
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FIG. 1: (Color online) Evolution of the three- and four-boson binding energies with the cutoff 
parameter A. 

The appearance/disappearance of trimer and tetramer excited states can be pointed- 
out in the A-dependence of the atom-dimer and atom-trimer scattering lengths a\2 and a\s 
shown in Fig. [2j when a new trimer (tetramer) state appears, ayi (013) has a discontinuity 
going to —00 and returning from +00. Thus, Figure [2] demonstrates that new trimer states 
appear around a§h = 28.12 and 31.7, whereas an additional tetramer state exists only for 
26.75 < aoA < 28.14. Much like the binding energies, the scattering lengths calculated with 
the rank 2 potentials show only weak A-dependence at aoA < 26.0, whereas more significant 
but monotonic A-dependence is seen in the rank 1 case. 

The EFT predicts universal nearly linear correlations between the three and four-particle 
binding energies revealed by the well-known Tjon line in the few-nucleon physics. As 
example we show in Fig. [3] the relation between trimer and tetramer ground state binding 
energies. Results obtained with n r — 1 potentials are roughly consistent with the findings 
of EFT in the whole considered A region where B± /B2 as a function of / B2 has slope 
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FIG. 2: (Color online) Evolution of the atom-dimer and atom-trimer scattering length with the 
cutoff parameter A. 

around 5. However, our n r = 2 results obviously violate linear correlations around ao-A = 
27.5 where the three-boson binding energies JBg are almost A- independent, whereas the 
four-boson binding energies Bf ) exhibit a very rapid increase. There is even a very narrow 
region around ciqA = 26.0 where trimer binding energies decrease, while at the same time 
tetramer binding energies grow with A. Nevertheless, nearly linear correlations are preserved 
at aoA < 26.0, also for the excited states. The linear correlation but with a different slope, 
around 2.4, takes place also at aoA > 28.5 where both B^ and B^ increase rapidly; 
however, due to small radii and high average kinetic energies, as will be shown later, this is 
not the regime where one should expect validity of the universal relations. 

Correlations between trimers ground and excited state binding energies are shown in 
Fig. |H In contrast to the nontrivial behavior of the results obtained with rank 2 potentials, 
contact interaction EFT, which ensures validity of the Efimov condition, proclaims univer- 
sality of the B i ™ ) /B 2 (B^ +1) /B 2 ) curve [3|]. The results obtained with realistic 4 He potentials 
lie slightly below the universal curve much like our rank 1 results do. On the other hand, 
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FIG. 3: (Color online) Correlation between trimer and tetramer ground state binding energies. 

the relation between B^ and B^ obtained with n r = 2 potentials agrees well with the rank 
1 curve only at g^A < 26.0 where binding energies decrease by increasing A. At a A = 26.0 
the .Bg /B 2 (B^ / B2) function exhibits a critical point behavior and recedes from the rank 
1 curve. The second excited trimer state appears at aoA = 28.12; the corresponding corre- 
lation function B§ / B 2 (B% j B 2 ) increases monotonically and is consistent with the rank 1 
S3 / B 2 (B% I B 2 ) curve till a^A = 29.2 but then deviates from it as well. 

Given such a nontrivial behavior of few-boson binding energies it is worth to study other 
bound state properties such as the expectation value of the kinetic energy operator, and 
the mean square radius, Rjy ■ As shown in Fig. [5j there is at least a qualitative correlation 
between kinetic and respective binding energies of the trimer and tetramer ground state, 
i.e., except for the vicinity of the critical points, Kffl increase with B^\ although the slope 
varies quite strongly. Evolution of the first excited state kinetic energies is more peculiar. 
In the region a^A > 26.0, where binding energies start to grow rapidly with A, the phase of 
the fast increase of kinetic energies is followed by the phase of fast decrease after which 
the K$ increases in a way similar to K^. 

The correlations between mean square radii R$ and the corresponding separation en- 
ergies [B^ — i%_J are shown in Fig. [6j We note that for the ground states < r = 
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FIG. 4: (Color online) Correlations of the trimer binding energies. 
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FIG. 5: (Color online) Correlations between the expectation values of kinetic energy and the 
corresponding binding energies. 
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FIG. 6: (Color online) Correlations between the mean square radii and the corresponding separation 
energies. 

0.144 do, thus, the observed deviations from the universal behavior are not surprising. As 
it is natural to expect, the radii decrease with increasing separation energy, often the de- 
pendence is not far from linear. The exception is the first excited state of the tetramer at 
aoA > 26.75. The origin of this peculiar behavior is unveiled in the next paragraph, when 
discussing the shape of non-local potential. Worth noting is that in a small window around 
aoA = 26.75, corresponding to the first rapid increase of B± , the first peak and the 
first i?4 minimum, the tetramers excited state kinetic energy is larger and its mean square 
radius is smaller that the corresponding quantities of the ground state, i.e., K± > K± and 

p(l) ^ p(0) 

Explanation of such a nontrivial behavior of the few-boson system is rather simple. For 
this aim one should simply take a look at the A-evolution of the nonlocal interaction, repre- 
sented in the configuration space, see Fig. [7J Potentials with small cutoff A have single wide 
and rather flat attractive plateau in the (r, r') plane. For larger A this plateau splits up into 
a rather complicated saddle-like surface with two asymmetric attractive regions separated by 
two symmetric repulsive regions. One narrow attractive region that deepens with increasing 
A is at the origin while wider and natter region is more distant. Due to large kinetic energies 
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required to squeeze particles into the narrow region at the origin, it is energetically prefer- 
able for the multi-boson system to reside in the distant attractive plateau. The resulting 
spectrum obeys the standard universality relations p, |6j, e.g., it has two tetramer states, 
one tightly bound and one weakly bound, which is essentially an atom weakly bound to the 
trimer in its ground state. However, at some point close to a A = 26.0 the attractive region 
at the origin becomes deep enough to accommodate the four-boson system, and some kind 
of a phase transition occurs. First, the excited state of the tetramer accommodates impor- 
tant part of its wave function closer to the origin thereby loosing its atom-trimer structure 
and experiencing a drastic reduction of the radius R± and a growth of the potential, ki- 
netic and binding energies. Further deepening of the small (r, r') region brings the tetramer 
ground state from the distant plateau to the attractive region close to the origin. The first 
excited tetramer state is pushed away to "recover the former place" of the ground state, 
thereby decreasing and increasing w± , and a second excited state appears in the place 
of the former excited state. As a consequence, at aoA = 27.0 the tetramer ground state is 
tightly bound in the narrow deep region close to the origin while excited states reside in the 
wider flatter region, the first one being tightly bound and the second one having an atom- 
trimer structure; furthermore, their relation to the trimer is nearly universal as for the two 
tetramer states in the region a^A < 26.0. Thus, at aoA = 27.0 it is the ground state of the 
tetramer with the radius ~ r /5 much smaller than the range of the interaction that is 
anomalous/non- universal one. Trimer has only three pair interactions (compared to six for 
the tetramer), therefore the corresponding phase transition in the three-boson system takes 
place at larger A values, when the attractive region near the origin is even more deepened. 
As the trimer ground state is brought to the narrow deep region as well, the tetramer states 
havo to change correspondingly. Since there is already tightly bound tetramer in that re- 
gion, the second excited tetramer state disappears and the first one acquires the atom-trimer 
structure. As a consequence, the separation energy B± — B% decreases while individual 
binding energies continue to increase and radii decrease, as can be seen from the last part 
of the i?4 curve in Fig. El 



Furthermore, we would like to mention that Delfino et al. 



20] have found a strong collapse 



of the trimer for a rank 1 separable potential with two-term form factors. Effect seen in that 



work may be well related with the one we discuss, since the potential of Ref. 



20j may be 



rewritten as a rank 2 potential with only two independent parameters and therefore has 
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the structure comparable to ours with several attractive and repulsive regions. However, 



in Ref. 



201 ] single two-body parameter, namely was fixed and thus collapse could be 



interpreted as a standard Thomas collapse, being accompanied by the rapid increase in 
the average kinetic energy of the dimer (or reduction of the effective range). On the other 
hand, our potentials are phase equivalent in a broad energy window, up to 50 B 2 , thereby 
constraining in addition two-body scattering length, effective range, and even higher terms 
of the effective-range expansion. Nevertheless, presence of finite effective range expansion 
terms (ro > etc.) does not prevent collapse in the multi-boson system. Indeed, our 
additional calculations indicate that varying upper limit of the energy included in the fit 
changes only the critical A values but not the qualitative picture of the collapse. 

One should probably question the existence of such "exotic" nonlocal potentials with two 
attractive regions. At present, we cannot give an example of a real physical system where a 
nontrivial behavior like described above is observed. However, we admit that some realistic 
nucleon-nucleon potentials H Q have even mOT e complex structure than that, containing 
several attractive regions. 

Finally, in Figs. El and [9] we show correlations between the atom-dimer (atom-trimer) 
scattering length and the trimer (tetramer) ground state binding energy. It is nearly linear 
for rank 1 results, reproducing the well-known Phillips line whereas rank 2 results agree with 
that line only up to aoA = 26.0 or even less in the case of a 13 . Such a behavior is consistent 
with previous findings of the present work. 



IV. CONCLUSIONS 



We calculated properties of three- and four-boson systems employing series of short-range 
two-boson potentials reproducing the same two-body binding energy and low-energy scat- 
tering observables. In a large part of the considered parameter space we obtain results that 
are consistent with the commonly established EFT results. However, a smooth variation 
of the potential parameters leads to the region with drastic deviations form the EFT pre- 
dictions, e.g., rapid variations of tetramer binding energies while trimer binding energies 
remain almost unchanged and the appearance of the second excited tetramer state, the rea- 
son for these phenomena being complicated structure of the potential with several attractive 
regions. 



13 




FIG. 7: (Color online) Configuration space (r, r') representation of the nonlocal rank 2 potentials 
in the form V = r(r\v\r')r' , obtained with A = 0.15 A -1 (top) and A = 0.27 A -1 (bottom). 
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FIG. 8: (Color online) Correlation between the atom-dimer scattering length and the trimer binding 
energy. 
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FIG. 9: (Color online) Correlation between the atom-trimer scattering length and the tetramer 
binding energy. 
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Lately models based on contact interaction EFT have found a huge success in explaining 
correlations between the few-particle binding energies in nuclear and cold-atom systems as 
well as establishing existence of universal relations between them. The results obtained by us 
demonstrate existence of strong deviations from the aforementioned universal behavior. In- 
deed, universal behavior is assured only if Efimov condition r <C R <C |ao| is satisfied, which 
is fulfilled if one employs contact-interactions (as is used to establish universal relations). 
Nevertheless, this condition is not guaranteed a-priori for multi-particle system. Moreover, 
the realistic N > 2 particle systems, even those which are believed to be ideally suited to 
apply EFT like cold atomic gases or nuclei, do not always satisfy the first Efimov condition 
r <C R. For example, r ~ R for 4 He trimer and tetramer ground states calculated with 
"realistic" potentials between 4 He atoms ll| or with the present toy model. The universal 
behavior still holds in case when the interaction is rather simple, having single attractive 
region and/or N > 2 and N = 2 systems are bound to the same attractive region, however, 
it is not guaranteed a-priori. A general N + 1 particle system may have very different states 
from those predicted by the universal EFT, even though all the spectra of N or less particles 
are well described. 

Presence of finite range (r > 0) does not guarantee that Thomas collapse would not occur 
in multi-boson system. Interaction can have short-range or off-shell structure that may not 
be seen with the low-energy probe. Such structure can have no or little effect for low-energy 
observables of the systems up to N particles, however may be exploited in the N + 1 boson 
system that, having more interacting pairs, can compensate larger kinetic energies and thus 
regroup itself into the shorter range domain. As discussed in Ref. 5] one can choose three- 
body scale for repulsive three-body interaction to regularize three-body system and prevent 
iV-particles from collapsing. Nevertheless, in special cases three-body low-energy observables 
may carry insufficient amount of information to set the physics beyond N = 3 system as we 
demonstrated in the present work and a four-body parameter may be needed 17| • However, 
this may take place only when the underlying potential has a nontrivial structure and the 
size of the few-body system R < rn becomes smaller than the range of the interaction. Thus, 

nn 

our results are complementary to the findings of Refs. |5|, |6J that four-body parameter is not 
needed under Efimov condition R» r assuring the universality. 
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